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We present new results for the amplitude A 2 for a kaon to decay into two 
pious with isospin 1 = 2: ReA 2 = l-50(4)stat(14)syst x 10“® GeV; ImA 2 = 

— 6.99(20)stat(84)syst X 10“^^ GeV. These results were obtained from two ensembles 
generated at physical quark masses (in the isospin limit) with inverse lattice spacings 
a~^ = 1.728(4) GeV and 2.358(7) GeV. We are therefore able to perform a continuum 
extrapolation and hence largely to remove the dominant systematic uncertainty from 
our earlier results [1, 2], that due to lattice artifacts. The only previous lattice com¬ 
putation of A —)• vrTT decays at physical kinematics was performed using an ensemble 
at a single, rather coarse, value of the lattice spacing [a“^ ~ 1.37(1) GeV]. We con¬ 
firm the observation reported in [3] that there is a significant cancellation between 
the two dominant contributions to ReA 2 which we suggest is an important ingredi¬ 
ent in understanding the AI = 1/2 rule, ReAo/ReA 2 — 22.5, where the subscript 
denotes the total isospin of the two-pion final state. Our result for A 2 implies that 
the electroweak penguin contribution to e'/e is Re(e'/e)EWP = —(6.6 ± 1.0) x 10“^. 
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I. INTRODUCTION 

Nonleptonic K —)■ tttt decays continue to be an important class of processes in the phe¬ 
nomenology of the standard model of particle physics. Historically it was in these decays 
that both direct and indirect CP-violation were discovered and the challenges for theoretical 
physicists include an explanation of the long-standing puzzle of the AJ = 1/2 rule and an ab 
initio computation of e'/e. Developments in the theoretical framework of lattice QCD and in 
efficient algorithms, together with the availability of the latest computing power, have made 
meeting these challenges feasible. A signihcant element of the current joint research program 
of the RBC and UKQCD collaborations is the evaluation of the K —)■ tttt amplitudes Aq and 
A 2 , where the subscript represents the isospin of the two-pion hnal state (which by Bose 
symmetry is restricted to 0 or 2). In this paper we present our latest results for A 2 . 

In [1, 2] we reported on the hrst results from a lattice determination of the amplitude A 2 
for K —)■ (7r7r)/=2 decays, where / is the total isospin of the two-pion hnal state: 

ReA 2 = 1.381(46)stat(258),yst 10“® GeV, Im^ = -6.54(46)stat(120)syst GeV . (1) 

This was the hrst quantitative calculation of an amplitude for a realistic hadronic weak 
decay and hence extended the framework of lattice simulations into the important domain 
of nonleptonic weak decays. As explained in the Introduction of [2], in order to obtain the 
result in Eq. (1) it was necessary to overcome a number of theoretical problems and exploit 
recent improvements in algorithms and the opportunities provided by increases in computing 
resources. The systematic errors in (1) are dominated by the fact that the calculation was 
performed at a single, rather coarse, value of the lattice spacing (a ~ 0.14 fm). We estimated 
these errors to be 0(15%). 

In this paper we repeat the calculation at two hner values of the lattice spacing and 
perform the continuum extrapolation.The simulations are carried out at physical pion masses 
(with unitary sea- and valence-quark masses) using our two new ensembles with lattice 
spacings a = 0.011 fm and a = 0.084fm. Our new result is presented in Eq. (63) and we 
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reproduce it here for the reader’s convenience: 


Re(A 2 ) = 1.50(4),tat(14)syst X 10-" GeV; Im(A 2 ) = -6.99(20)stat(84)syst x IQ-^" GeV. 

( 2 ) 

A very interesting feature of our earlier calculation of A 2 was the observation that the two 
dominant contributions to ReA 2 show a signihcant numerical cancellation [3]. We argued 
in [3] that this cancellation is an important element in the explanation of the AI = 1/2 
rule, ReAo/ReA 2 ~ 22.5. We conhrm this cancellation in the present calculation. Of course, 
before we can claim that we fully understand the A/=l/2 rule, we need to compute Aq 
at physical quark masses and momenta; this calculation is even more challenging than the 
evaluation of A 2 but is under way. For the status of this calculation we refer the reader 
to [4]. 

The structure of the remainder of this paper is as follows. In the next section we present 
the parameters of the two ensembles used in this calculation. The evaluation of the bare 
matrix elements and the renormalization of the lattice operators are discussed in Secs. Ill 
and IV respectively. We consider hnite-volume effects in Sec. V and present an overview 
of the different sources of systematic uncertainty in Sec. VI. We perform the continuum 
extrapolation in Sec. Vll and present our hnal result in Eq. (63). Section Vlll contains 
our conclusions and a brief discussion of the prospects for the reduction of the errors in 
A2 as well as for the calculation of Aq. There is one appendix in which we reproduce the 
calculation from [5] of the Lellouch-Liischer factor for hnite-volume corrections in the context 
of chiral perturbation theory. This calculation demonstrates how to disentangle the hnite- 
volume corrections which decrease exponentially with increasing lattice volume (a source of 
systematic error) from those which decrease as a power of the volume (which are corrected 
by the Lellouch-Liischer factor). This calculation also clarihes a misunderstanding of these 
effects in the literature [6]. 

II. DETAILS OF THE SIMULATION 

The calculations described below have been performed on two new 2-1-1 havor ensembles 
generated with the Iwasaki gauge action and with Mobius domain-wall fermions [7] . The 
parameters of the ensembles are 
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(i) 48=^ X 96 X 24 with /3 = 2.13 = 1.728(4) GeV); 

(ii) 64^ X 128 x 12 with (3 = 2.25 (a-^ = 2.357(7) GeV). 

These two ensembles use the Mobius variant of domain wall fermions [8] with a Mobius 
scale factor a = 2. For compactness of notation we will refer to these ensembles as 48^ and 
64^ respectively. The lattice spacing and quark masses were set by choosing the masses of 
the pion, kaon and the G-baryon to be equal to their physical values. The corresponding 
sea-quark masses are arriud = 7.8 x 10“"^ and arris = 3.62 x 10“^, with the residual mass 
anires = 6.19(6) x 10“^ for the 48^ ensemble and arriud = 6.78 x 10“^, arris = 2.661 x 10“^ 
and arrires = 2.93(8) x 10“^ for the 64^ ensemble. The two ensembles have approximately the 
same physical volume with spatial extent L ~ 5.5 fm, enabling the continuum extrapolation 
to be separated from hnite-volume effects which we estimate separately. For more details on 
these ensembles see [7] and we will return briefly to the determination of the lattice spacings 
in the context of the continuum extrapolation in Sec. VII. 

The results presented below were obtained using 76 gauge conhgurations on the 48^ en¬ 
semble and 40 on the 64^ ensemble. The large statistical uncertainty one expects with a 
relatively small number of gauge conhgurations can be signihcantly reduced if we perform 
many measurements on each conhguration in which the sources and sinks are simply trans¬ 
lated in space and time [7]. Performing multiple measurements on the same conhguration 
ohers two important opportunities for increased efficiency. First if we can use a low-mode 
dehation method such as eigGG [9] we will be able to amortize the setup costs of such an 
approach over a large number of inversions. Second we can use the all mode averaging 
technique [10] and perform most of these many inversions at reduced precision and use a 
relatively few accurate inversions to determine a correction that guarantees systematic dou¬ 
ble precision but with an additional (usually small) statistical error that rehects the small 
number of accurate solves. Specihcally for the 48^ ensemble, the eigGG method was used 
in single precision with 600 approximate low-lying eigenvectors and a stopping residual of 
10“^. The approximate (wall source) propagators were computed on all 96 time slices. The 
accurate solves used to correct the approximation were computed on time slices 0, 76, 72, 
68, 64, 60 and 56 with Gonjugate Gradient (GG) stopping residual 10“®. (This choice of 
time-slice separations is not related to the 77 —)■ tttt calculation presented here but to an 
accompanying calculation of Bx [7].) To ensure that no bias results from the choice of 
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rriK 

E 

-^TTTT 

TflK -^TTTT 

48^ (lattice units) 
64^ (lattice units) 

8.050(13) X 10-2 
5.904(14) X 10-2 

2.8867(15) X 10-^ 
2.1531(14) X 10-1 

2.873(13) X 10-1 
2.1512(68) X 10-1 

1.4(14) X 10-3 
9(10) X 10-^ 

48^ (MeV) 

643 (MeV) 

139.1(2) 

139.2(3) 

498.82(26) 

507.4(4) 

496.5(16) 

507.0(16) 

2.4(24) 

2.1(26) 


TABLE I; Pion and kaon masses and the 1=2 two-pion energies in lattice and physical 
units measured on the 48^ and 64^ ensembles. The momentum of each of the hnal-state 
pions is ±7r/L in each of the three spatial directions. 


inexact solves for which the correction is calculated, this complete pattern of source time 
slices for the accurate solves was shifted by a different random time displacement on each 
conhguration. A similar procedure was used on the 64^ ensemble but with 1500 low modes 
and a stopping residual of 10“^ for the approximate solves and accurate solves on time slices 
0, 103, 98, 93, 88, 83, 78 and 73. On both ensembles, the accurate CG solves were also com¬ 
puted using eigCG, exploiting the approximate eigenvectors created during the inaccurate 
applications of eigGG. 

Measurements on the 48^ and 64^ ensembles are separated by 20 and 40 molecular dy¬ 
namics (MD) units respectively. In order to study the effects of autocorrelations we bin the 
data. We hnd that the effects are small, typically leading to a variation of the statistical 
errors of less than 10%. The results presented below were obtained after binning the 76 
conhgurations of the 48^ ensemble into 19 bins of 4 conhgurations and the 40 conhgurations 
of the 64^ ensemble into 8 bins of 5 conhgurations. The 40 conhgurations from the 64^ en¬ 
semble are precisely those used in the global analysis reported in [7]. The 76 conhgurations 
from the 48^ ensemble include 73 of the 80 used in [7]. We have however, repeated the 
relevant analysis of [7], including the determination of the lattice spacings, using precisely 
the 76 conhgurations for which we have computed A 2 . This makes it possible to compute 
standard jackknife errors for our physical results which necessarily depend upon the value 
of the lattice spacing. 

The pion (m^) and kaon masses (mx) as well as the energies of the 1 = 2 two-pion state 
obtained on the two ensembles are shown in Table I. The htting ranges used for pion 
and kaon masses as well as two pion energies were from 10 to 86 on the 48^ ensemble and 
from 10 to 118 on the 64^ ensemble. These choices were motivated by the plateaus in the 
ehective mass plots shown in Figs. 1-2. The ehective mass of the kaon, is dehned 
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FIG. 1: Effective mass plots for the kaon correlation functions on the 48^ ensemble (left) 

and 64^ ensemble (right). 
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FIG. 2: Effective mass plots for the two-pion correlation functions on the 48^ ensemble 

(left) and 64^ ensemble (right). 


numerically by the ratio 


Cxit + l) cosh(m|(?'(t + 1 — T/2)) 

Cxit) cosh(m|f(t — T/2)) 

and the two-pion effective mass, E®®, is found by inverting 

+ 2) — G,r7r(t -l- 1) -f- 

+ 1 ) - ^ - 1 - 


The two-point correlation functions Ck and Ct^tz are dehned explicitly in Eq. (22) below 
and the differences in the numerator and denominator on the left-hand side of Eq. (4) are 
introduced to eliminate the constant C in Eq. (23). 
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The pion and kaon masses correspond closely to their physical values. We will explain 
below that the pions are given a momentum tt/L in each of the three spatial directions and 
from the table we see that with this choice ~ mx and the K —)■ tttt matrix elements 
correspond to the on-shell (within statistical errors) decay of a kaon in the center-of-mass 
frame. We now discuss the evaluation of the matrix elements. 


III. EVALUATION OF THE BARE MATRIX ELEMENTS 

A —)■ TTTT decay amplitudes are defined by 

= ((7r7r)7=2,o I Hw \ A°), (5) 


where Hw is the component of the weak Hamiltonian which changes the strangeness by one 
unit. The weak Hamiltonian can be separated into short and long distance contributions by 
using the operator product expansion: 

Hw = 5 ^ cApmp), ( 6 ) 

^ i 

where Gp is the Fermi constant, Ws and Vud are Cabibbo-Kobayashi-Maskawa (CKM) ma¬ 
trix elements, the Qi are all the possible dimension-6 operators which contribute to the 
decay and Ci are the corresponding Wilson coefficients which contain information about the 
short distance physics. The Ci take the form Ci = Zi + ryi where r is the ratio of CKM 
matrix coefficients r = — . 7%)'^ . 

V*sVud 

In this paper we only consider AJ = 3/2 decays where the two-pion hnal state has total 
isospin 2. The nonperturbative contribution to the decay amplitude is contained in the 
matrix elements: 


- 


= I I I Q2Y= 


K+ _ 


.AI=3/2 

'Al3=l/2,i 


There are only three operators which contribute to A 2 , which we label according to their 
chiral SU{3)l x SU{3)r transformation properties. We have one (27,1) operator and two 
electroweak penguin operators labeled (8,8) and (8, 8)mx, where the subscript mx denotes a 



color mixed operator. Explicitly, the operators are given by 


^(27 1 ) ( S ) 

^(8,8) ~ {sidi)L {jdjUj — djdj ^+ (siUi)liujdj)n, (9) 

^(8,8)mx i,^idj')L (ujU-i djd^^ -\- (ySj^Uj)ij{Ujdj^ji. (^0) 

The subscripts L and R denote the left- and right-handed spin structures respectively: 

(gig 2 )L = - 7^)92 and (gig 2 )R = gi7^(l + 7^)92- (11) 

The Lorentz indices are understood to be contracted between the two parentheses in each 
of the operators in Eqs. (8)- (10) and f, j are color indices which are summed from 1 to 3. 

Below we will conhrm the feature found in our earlier work[l, 2] that the dominant 
contribution to Re(R 2 ) comes from the (27,1) operator, while the dominant contribution to 
Im(R 2 ) in the MS scheme at 3 GeV comes from the (8, 8)mx operator. We can now write the 
expressions for the A 2 amplitude, which are 



The relative factor between the two expressions is due to the different Clebsch-Gordan 
coefficients. 

A major challenge in the calculation of A 2 (and even more so in the calculation of Aq) is 
to ensure that the pions have physical momenta. In the center-of-mass frame with periodic 
boundary conditions, the ground state for the two-pion system has each pion at rest. The 
evaluation of matrix elements at physical kinematics therefore corresponds to the contribu¬ 
tion from an excited two-pion state resulting in a considerable loss of precision. We can 
avoid the necessity of multiexponential hts to extract the excited state contribution by uti¬ 
lizing the technique suggested in [11, 12] and applied successfully in our original calculation 
of A2 [1, 2]: we introduce antiperiodic boundary conditions for the (valence) d-quark in all 
three spatial directions, and periodic boundary conditions for the u- and s-quarks [11]. We 
then exploit the Wigner-Eckart theorem to relate matrix elements to those for 
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the unphysical transition —)■ . The relation is 

flTTTry, I QaSv2 ,» I \ QlHS/2,i I K+) • (13) 

-^((7r+7r0| + (7r07r+|) (tt+tt+I 

The indices I and Is label the two-pion state’s total and third component of isospin respec¬ 
tively. With antiperiodic boundary conditions in three spatial directions, the ground 

state has total momentum 0, with each pion having momentum |pV| = It can be 

seen from Table I that is very close to rriK on both the 64^ and 48^ ensembles. (For the 
smaller physical volume in our original calculation [1, 2], we imposed antiperiodic boundary 
conditions for the d-quark in two spatial directions in order to achieve ~ ttik-) Note 
that with both periodic and antiperiodic boundary conditions on the d-quark, the lowest 
momentum of the 7r° meson is zero; this motivates the use of the Wigner-Eckart theorem to 
reformulate the calculation to that of a matrix element with a hnal state. 

The operators Qa/;((^ 3/2 '''’hich appear on the right-hand side of Eq. (13), and which cor¬ 
respond to the Q^j~^y 2 operators in Eqs. (8)- (10), are 

Q(‘21,1) Lij^jdj') Li Q(8,8) ('Sjdj) 2 y('Ujdj)/j, Q(8,8)mx ('^idj')^ (lijdj)^. (14) 

To simplify the notation we have dropped the labels AJ = 3/2 and A/^ = 3/2 on the 
operators in Eq. (14); this will be implicit in the following. In this paper we compute the 
K —)■ TTTT matrix elements of the three operators in Eq. (14). 

The factor of 3/2 in Eq. (13) is a combination of \/3/2 coming from the Clebsch-Gordan 
coefficients and the Wigner-Eckart theorem, and a further \/3 corresponding to the simple 
choice for the normalization of operators in Eq. (14). The amplitude A 2 is given in terms of 
the matrix elements Mi by 

.42 = ^ Z 0(f‘) V.. (15) 

Since it is the —)■ matrix elements which we compute directly in this paper, we 

choose the compact notation Mi = . The label i runs over the three operators in 

Eq. (14). 



10 


A. Evaluation of the correlation functions 

The bare matrix elements are obtained from the computation of two- and three-point 
correlation functions. The three-point functions are 

= (0 I Q^{top) al,(0) | 0), (16) 

where Qi is one of the three operators in Eq. (14) and ax and are interpolating operators 
for the kaon and two-pion state respectively. For ax and we take Coulomb gauge-hxed 
wall-source operators dehned as follows: 

^K{t) = ^ s(fi,t) 7®M(T2,t), (17) 

Xl,X2 

= [d{t)'y^u{t)] [d{t)'y^u{t)] , (18) 

where in (18) we have used the cosine momentum sources for the d-quark: 

d(t) = ^ d{x, y, z, t) cos{xp^) cos{ypy) cos{zp;.). (19) 

x,y,z 

d{x, y, z, t) represents the d-quark held and the components of momenta satisfy Px = Py = 
Pz = TijL. Just as for the w-quark source in Eq. (17), the w-quark sources in a-x-,, shown in 
Eq. (18) are given zero momentum by summing them over the full spatial volume, evaluated 
in the Coulomb gauge. As explained in Ref. [2] the cosine source described above creates 
d-quarks with both signs for each component of the three momentum ±pj, for i = x, y and 
2 ;. This will then produce pairs of pions with total momentum in each direction of ±27r/L 
in addition to the desired value of 0. For the three-point functions described in Eq. (16), 
the zero total momentum of the decaying kaon and three-momentum conservation imply 
that the nonzero tt-tt momenta cannot occur. For the two-point function dehned in Eq. (22) 
below we use a tt-tt sink which is different from the source and which explicitly projects 
onto TT-TT states with zero total momentum, as described in Ref. [2]. A further subtlety, 
not described in that reference, relates to the possible angular momentum of the two-pion 
state. For our two identical tt’*' bosons which carry equal but opposite momenta, there 
are actually four possible states given our boundary conditions. Specihcally, the tt’*' which 
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FIG. 3: Diagrammatic representation of the K —)■ tttt three-point fnnction dehned in 
Eq. (16). The strange-qnark propagator is explicitly labeled, the remaining lines represent 

light-qnark propagators. 


carries px = +t^/L may have four possible values for the other momentum components: 
Py = ±71/L and Pz = ivr/L. These four states form a four-dimensional representation of the 
cubic symmetry group, which decomposes into two irreducible representations: a singlet (Ai) 
and a triplet (T 2 ), out of which only Ai contains an s-wave contribution. Since the lowest 
energy level of the finite-volume 1 = 2 s-wave spectrum of the Ai representation is nearly 
degenerate with the lowest energy level of the d-wave spectrum of the T 2 representation, it is 
important that we use the cubically symmetrical source specified in Eq. (19) which couples 
only to the Ai state of interest. 

The spinor and color labels are contracted within each set of square parentheses in 
Eq. (18). A schematic diagram of the correlation function is shown in Fig. 3. 

We have evaluated C^^'^'^{top) for a range of values of the source-sink separations 
For the 48^ (64^) ensemble we performed the calculations for values of between 24 and 
39 (26 and 36). These separations were chosen to be large enough for the plateau region 
to give a reliable fit and small enough for the around-the-world effects to be small. The 
fitting ranges were chosen to be from 10 to — 10 for both ensembles. These choices are 
motivated by the locations of plateau regions in Fig. 4. 

For sufficiently large time separations Gp and t,r 7 r — Gp, the expected time dependence of 
C'f^’^^(top) is 


( 20 ) 
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FIG. 4: iF —)■ TTTT three-point correlation function on the 48^ lattice (left) and 64^ lattice 
(right) with a kaon-pion separation of = 26. 


where 

= K^TTr I (T^^(O) I 0)1 and iVi^-= |(iF| 0 -^( 0 ) |0)|. (21) 

We have introduced the label “bare” as a reminder that are matrix elements of the 

bare operators in the lattice regularization which we are using. The renormalization of the 
operators is discussed in the following section. For illustration, in Fig. 4 we plot 
computed on each of the two ensembles for = 26. The observed plateaus are a manifes¬ 
tation of the fact that the volumes have been tuned so that — 'mx [cf. Eq. (20)]. 

We obtain the matrix elements Mj by htting Eq. (20), using the values of Nk, ttik 

and obtained from htting (under the jackknife) the correlation functions, 

G^^(f) = (0 \al^{t,p = 0) a^^(0)| 0) and Cxit) = (0|ax(t)cxx(0)|0), (22) 

which have the following time dependence: 

+ C) , (23) 

t^oo ' 

Cxit) -^ \Nk\^ . (24) 

i —>-00 ' ' 

The “f oo" limit should be understood as taking a sufficiently large time separation so 
that excited state contributions are negligible. Introducing the constant C in Eq. (23) allows 
one to account for possible around-the-world effects in C-,,.,,. 
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FIG. 5: Ratios of iF —)■ vrTr three-point correlation function to the two-point functions 
(Eq. (25)) on the 48^ lattice (left) and the 64^ lattice (right) with a kaon-pion separations 

of fjrvr = 27 and 36 respectively. 




3 /Cfbare 


48^ ensemble 
64^ ensemble 

3.700(35) X 10-^ 
1.371(11) X 10-^ 

9.171(69) X 10-3 
3.942(39) X 10-3 

3.058(23) X 10-2 
1.308(13) X 10-2 


TABLE II: Results for the bare —)• vr’^Tr^ matrix elements in lattice units. Only 

statistical errors are shown. 

As a check, we can also construct the time-independent ratio of the correlation functions: 

_ MtL^ f25) 

Cxit) C'jr7r(G7r “ t) 

This ratio is plotted for tTrn = 26 in Fig. 5. As anticipated, all three operators exhibit a 
constant behavior in the region where the contribution from excited states is negligible. 
Equation (25) is expected to hold in the region 0 Gtt ^ T, where T is the total time 

extent of the lattice. In this region “around-the-world” effects arising from different time 
orderings of the operators can be neglected. 

The values of the bare K~^ —)■ matrix elements are shown in Table II. The entries 

have been obtained by performing weighted averages (under the jackknife) over the values 
obtained for each choice of 
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IV. RENORMALIZATION OF THE OPERATORS 

Having determined the matrix elements of the bare operators in the lattice regularization 
we now have to combine them with the remaining factors in Eq. (6) to obtain A 2 . The Wil¬ 
son coefficients [Ci(/r)] and composite operators appearing in Eq. (6) are separately 

renormalization scheme and scale (/i) dependent. To obtain the physical amplitudes they 
must be combined in the same scheme and at the same scale. The Ci{n) are calculated in 
perturbation theory for which it is convenient to use the MS -NDR scheme (called MS in 
the following). NDR stands for “naive dimensional regularization” prescription for the 7 ® 
matrix, which preserves the anticommutation relations with other gamma matrices [30]. The 
matrix elements calculated in Sec. Ill, on the other hand, were obtained using bare operators 
with the lattice spacing as the ultraviolet regulator with the lattice discretization of QCD. 
The operators can be renormalized nonperturbatively, but only into schemes for which the 
renormalization condition can be imposed on lattice Green’s functions. The MS scheme, 
which is based on dimensional regularization cannot be simulated in a lattice computation. 
Our procedure is to start by renormalizing the operators non-perturbatively into schemes 
which can be simulated, specihcally the “regularization-independent symmetric momentum” 
(RI-SMOM) schemes [13] as described in detail in [2] and briefly summarized below. The 
matching between the RI-SMOM and MS schemes is necessarily performed in perturbation 
theory and is currently known at one-loop order. (Below we also present the matrix elements 
in two RI-SMOM schemes so that if the perturbative coefficients are calculated to higher 
order in the future, these matrix elements can be used to reduce the systematic uncertainty 
in A 2 due to the truncation of the perturbation series.) 

We now briefly summarize the renormalization procedure. We write the hve-point am¬ 
putated Green’s functions of the three operators in Eq. (14) as a three-component vector 
A = (A( 27 q), A( 8 , 8 ), A( 8 ^ 8 )jnx) = (Ai,A 2 ,A 3 ), and impose a renormalization condition of the 
form 

P{A 7 ,.)}=F, (26) 

where P is a vector of projectors and F the corresponding tree-level matrix. Denoting the 
tree-level contribution by the superscript ( 0 ) and including explicitly the spinor and color 
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labels, the matrix F is given by 


F 



rp,BA-,DC 
I (3a.\5') 



AB]CD 

q;/3;7(5 


tj- 


(27) 


Here greek letters label spinor components, the uppercase roman letters represent color 
indices and i,j = 1,2,3 denote the operators and projectors. For illustration, the tree-level 
value of the Green’s function of is 


A 


( 0 ) 

1 


AB-,CD 

a/3;'yS 


[WafiWyS + ( 7 ^ 75 )a/3 (7''75)75] 

- + (7^75)a5(7^75)7/3] ■ 


(28) 


For the renormalization we only consider the parity-even component of the four-quark op¬ 
erators. 

The choice of projectors is not unique and we implement two different sets known as the 
7 ^ and ^-projectors, given explicitly by 




JI-,LK 

I3(y;5'y 


^ + (7V)/3a(7V)57] ^ 

[(7^)/3a(7'')57 - (7''7^)/3a(7V)57] 


and 


' pd)' 


JI;LK 

^a^;S'y 


M)M)5'y - {iA)pa{iA)5^] 
y M0^{4)sa - {iA)0't{iA)sa] y 


The corresponding matrices F read 


p(7^) 


728iV(]V + l) 0 o' 
0 12SN^ 12SN 

^ 0 12SN 12SN^ I 


(29) 


(30) 


(31) 
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and 



^ 32N{N + 1) 
0 

V 0 


0 0 ^ 
32N^ 32N 
32N 32N‘^ j 


(32) 


where iV = 3 is the number of colors. 

The hnal result for the amplitude is, of course, independent of the choice of intermediate 
scheme dehned by P, but comparing the results obtained with different projection operators 
gives us an estimate of the systematic uncertainty due to the truncation of perturbation 
theory in relating the RI-SMOM schemes to the MS schemes. 

The renormalized operators are related to the bare ones by a matrix relation of the form 


Qfip) ^ Z„{pa)Q]'-’(a). 


(33) 


In order to extract the renormalization constants we follow the standard procedure [14, 15] 
and compute numerically the amputated Green’s functions of the bare operators in Eq. (14) 
with particular choices of external momenta (as discussed below) on Landau gauge-hxed 
conhgurations. We next solve Eq. (26) which we rewrite in the form 

= (34) 


where is the quark field renormalization constant and /i is the renormalization scale, 

which we ultimately choose to be 3 GeV. 

The choice of Zq is also not unique, and we use the following two cases: 



Tv 


tvao^, 


and 



1 

48 


1VAS)7", 


(35) 


where A.y is the three-point amputated Green’s function of the local vector current and Zy 
is the renormalization constant of the local vector current. In practice, we multiply each 
side of Eq. (34) by the square of the corresponding side of Eq. (35). This eliminates Zq and 
after this multiplication the left-hand side of Eq. (34) contains the ratio of renormalization 
factors ZijjZy. Zy is then calculated by imposing the Ward identity Zy{P | V | P) = 2mp, 
where is the local vector current and |P) is the state of a pseudoscalar meson P at rest 
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FIG. 6 : Momentum flow defining a renormalization condition of a four-quark operator in 
RI-SMOM scheme. The momenta are chosen so that pl= pl = (pi — ^ 2 )^ = 


with mass mp; this is explained in detail in [7]. 

The choice of projection operator for the four-quark operator and Zq dehnes a renormal¬ 
ization scheme, which we will label (a, h) with a, 6 G 7 ^^, ^ for the choice of and In 
particular, we consider the ( 7 ^, 7 ^^) and (^,^) schemes, having found in earlier studies that 
the perturbative conversion to the MS scheme is more precise in these schemes. This is based 
on the observation that the nonperturbative running is generally closer to the perturbative 
one for these schemes for the four-quark operators in Eq. (14) [ 2 , 16]. As explained below, 
we follow our previous practice and choose the (^,^) scheme for our central value and the 
( 7 ^, 7 ^) scheme to estimate the error due to the perturbative conversion to the MS scheme. 

Chiral symmetry suppresses mixing of operators in different irreducible representations of 
the chiral symmetry group, so that if the symmetry is exact, Zij is a block diagonal matrix 
with a 1X 1 block corresponding to the renormalization of the (27,1) operator and 2x2 block 
corresponding to the mixing of (8, 8) and (8, 8)mx operators. In a massless renormalization 
scheme with a chiral discretization such as the domain-wall action, we expect a mixing 
pattern very similar to this, but with a small 0{{amresY) mixing between the blocks. 

The mixing of the operator Q{ 27 ,i) with either of Q{s,s) or Q( 8 , 8 )mx due to explicit chiral 
symmetry breaking induced by hnite Lg is proportional to {am^esf' (which is ^ 3.6 x 10“^ 
in this work). Such mixing can result from two mechanisms [17, 18]. First, both quarks in 
a left-handed q-q pair in Q( 27 ,i) can propagate in the hfth dimension from the left-hand to 
the right-hand wall, exploiting numerous but exponentially damped modes which even in 
perturbation theory link the left- and right-hand walls. This will change the (27,1) operator 
into one transforming as the ( 8 , 8 ) representation, but requires the propagation of two quarks 
from the left-hand to the right-hand wall. This incurs a penalty of (arures)^ since one power 
of the residual mass results from the hfth-dimensional mixing of the left- and right-handed 
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components of a single quark. 

The second mechanism is nonperturbative and more subtle. For this case the propa¬ 
gation results from the left-right tunneling that can be caused by an eigenvector of the 
five-dimensional transfer matrix with a near-unit eigenvalue. Such eigenvectors permit 0(1) 
left-right mixing but are rare and therefore give a small contribution to rrires- Under some 
circumstances such modes can simultaneously allow a number of quark flavors to flip chi¬ 
rality. However, to change a (27,1) representation into an (8,8) one, both a quark and an 
antiquark must flip chirality which requires two distinct transfer matrix eigenvectors and is 
therefore also doubly suppressed by a factor {aniresY- Such doubled suppression will not 
occur for the mixing between the operator Q{ 27 ,i) and, for example, an operator in the (6, 6) 
representation. Here a single transfer matrix eigenvector with near-unit eigenvalue can re¬ 
sult in a 0(ami.es) mixing between Q( 27 ,i) and (s(l -|- 'y^)d) (m(1 -|- 7^)(i) by allowing both a 
u- and a d-quark (localized near this eigenvector) to flip chirality. This kind of mixing has 
been studied for example in [19] and it was found to be largely suppressed by our choice of 
kinematics, as explained below. 

In order to suppress physical infrared chiral-symmetry breaking effects we choose to 
impose the renormalization conditions with the kinematics indicated in Fig. 6 with pi = = 

{Pi ~ P 2 Y = We compute the Green’s functions for several momenta and interpolate to 
p = 3 GeV using a quadratic Ansatz. Using partially twisted boundary conditions, we have a 
good resolution around the targeted momentum. The momenta in such RI-SMOM schemes 
are chosen so that there are no “exceptional” channels, i.e. no channels in which the square 
of the momenta is small [13]. (This is in contrast with the original RI-MOM scheme [14, 15] 
in which pi = p 2 .) We have already checked that with do main-wall fermions and this 
choice of kinematics the chirally forbidden matrix elements are numerically negligible [2]. 
In the present computation, we use the 48^ and 64^ ensembles which have physical light 
and strange sea-quark masses. However, the light-quark mass is used in all of the valence- 
quark propagators in the hve-point Green’s functions, including those for both light and 
strange quarks. We do not extrapolate either the sea- or valence-quark masses to zero and, 
strictly speaking, do not work in the chiral limit. In practice the light-quark masses are 
sufficiently small that their effects are negligible as is the nonzero mass of the strange sea 
quark. Gomparing our results with those of our previous work (with Shamir domain-wall 
fermions) where a chiral extrapolation was performed we hnd agreement at the per-mille 
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level or better. 

We find that all the chirally forbidden renormalization factors are smaller than 10“^, so 
we set the corresponding matrix elements of Pi{Aj} to zero and finally obtain the renormal¬ 
ization matrices: 


= 3GeV) 


= 3GeV) 


^0.4617(3) 0 

0 0.5302(4) 

^ 0 -0.0386(1) 

^0.4822(3) 0 

0 0.5305(4) 

^ 0 -0.0637(1) 


0 ^ 

-0.07018(6) 
0.4451(5) ^ 

0 ^ 

-0.07135(7) 
0.5052(6) j 


for the 48^ ensembles and 


zj^ASiP = 3GeV) 


z'lfJLip = 3GeV) 


^0.5194(2) 

0 

^0.5399(2) 

0 

V 0 


0 

0.5774(2) 

-0.02797(7) 

0 

0.5782(2) 

-0.05230(4) 


0 ^ 

-0.0751(1) 
0.4431(6) j 

0 ^ 

-0.0761(1) 
0.4990(5) j 


(36) 


(37) 


(38) 


(39) 


for the 64^ ensembles. With momentum sources [20], only a few configurations are needed to 
obtain an excellent statistical precision. The number of Landau gauge-fixed configurations 
used to obtain these results varies between 5 and 15. The statistical errors were estimated 
with 200 bootstrap samples. The matrices in Eqs. (36) ~ (39) are the ones used in Eq. (33) 
to obtain the operators renormalized in the RI-SMOM schemes at the scale /i = 3 GeV from 
the corresponding lattice bare operators. 

The procedure described above enables us to calculate the matrix elements of the op¬ 
erators in Eq. (14) in the (continuum) RI-SMOM schemes with a very small systematic 
uncertainty due to the renormalization. The Wilson coefficients however, are computed in 
the MS scheme and so we have to match the RI-SMOM schemes to the MS one. We repeat 
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48^ ensembles 

64^ ensembles 

Re(A2) (7^7^) 
Im(A2) (7'',7'') 
Re(A2) (i^) 
Im(A2) (i^) 

1.346(ll)stat(l)NPR X 10“® GeV 
-5.739(46)stat(8)NPR X 
1.386(12),t,t(l)NPRX lO-^GeV 

—6.174(49)stat(9)NPR X 10“^^ GeV 

1.4029(93)stat(ll)NPR X 10-SGeV 
-6.143(73)stat(9)NPR X 10-13 GeV 
1.4386(95)stat(ll)NPRXl0-3GeV 
-6.548(78)3tat(10)NPR X IQ-i^ GeV 


TABLE III; The amplitude A 2 calculated using two different intermediate RI-SMOM 
schemes. The two errors, labeled “stat” and “NPR” are the statistical uncertainties in the 
evaluation of the bare matrix elements and respectively. Discrepancies in the results in 
the two schemes are attributed to the truncation in the matching to the MS scheme. 


that this matching is perturbative and at present is only known to one-loop order [21]; this 
limitation amplihes the uncertainty due to the renormalization. This uncertainty could be 
reduced by extending the perturbative calculations to higher orders. Future lattice calcula¬ 
tions could also help here by using step scaling to run the renormalization constants obtained 
in the RI-SMOM schemes nonperturbatively to larger momentum scales. The perturbative 
matching to the MS scheme can then be performed at these larger scales where the cou¬ 
pling constant is smaller, leading to smaller uncertainties. We now estimate the current 
uncertainty due to the matching. 

To estimate the uncertainty due to the truncation of the perturbative matching factors, 
we note that the matrix elements in the MS scheme should be independent of the choice of 
intermediate RI-SMOM scheme. Differences in the results are observed (see Table III) and 
attributed to the truncation. Following the procedure in [1, 2] we take the result obtained 
using the (^, intermediate scheme as our central value and the difference of the results 
obtained using the two schemes as an estimate of the systematic error. This uncertainty 
is marked as “NPR (perturbative)” in the error budgets presented in Tables IX and X in 
Sec. VI. The uncertainties marked as “NPR (nonperturbative)” are the statistical errors in 
the evaluation of Z^j. 


V. FINITE-VOLUME EFFECTS 


The presence of two pions in the final state in iP —)■ tttt decays leads to hnite-volume 
corrections which decrease as inverse powers of the volume, in addition to the exponential 
correction present in simpler quantities such as decay constants and form factors. The power 
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■^TTTT 

q 

6 (radians) 

as 

dq 

dip 

dq 

48^ 

0.2873(13) 

0.9087(61) 

-0.158(22) 

-0.174(24) 

3.7147(20) 

643 

0.21512(68) 

0.9157(43) 

-0.184(16) 

-0.201(17) 

3.7171(15) 


TABLE IV; Contributions to the Lellouch-Liischer factor on the 48^ and 64^ ensembles. 
The rate of change of the phase shift was calculated by using a linear approximation in 

momentum as explained in the text. 


corrections result in a multiplicative correction to the matrix element [22]: 


(tttt I Hw I K)oo = F{nn \ Hw \ K)fv ■ 


(40) 


The subscripts cxd and FV correspond to inhnite and hnite volume respectively, and the 
factor F is given by the Lellouch-Liischer formula [22]: 


= STiq + 

\oq oq 


86 \ ttikE, 


2 

TTTT 


po 


(41) 


where p is the magnitude of the momentum of a pion in the center-of-mass frame given by 
P = \/ ~ q is dehned as g = pL/27r. Since the mesons satisfy antiperiodic 

boundary conditions in all three spatial directions, the function 0 in this case is defined by 
the condition: 


tan0 


gTT^/^ 

^oo(i; q) ’ 


■^oo(i; q ) 


(42) 


6 is the two-pion s-wave phase shift, which can be calculated using the Liischer quantization 
condition, 5{q) -f 0(g) = rni, but the calculation of the derivative in Eq. (41) requires an 
approximation. 

The results presented in Table IV were obtained using the approximation that <5 is a linear 
function of the momentum between 0 and p. Since the second term in the parentheses on the 
right-hand side of Eq. (41) is much smaller than the first and given the remaining systematic 
uncertainties discussed in Sec. VI, this procedure gives an adequate approximation. In order 
to estimate the error due to this approximation we also evaluate the derivative using the 
phenomenological curve of Ref. [23] illustrated in Fig. 7; we take the difference of the two 
procedures as an estimate of the corresponding uncertainty. For our central value we use 
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FIG. 7: Comparison of / = 2 two-pion s-wave phase shifts calculated using Liischer’s 
formula with the phenomenological curve from Ref. [23]. The computed results are 
consistent with the phenomenological curve. 


the linear approximation for the derivative of the phase shift so that it is independent of 
phenomenological estimates. 

At the pion momentum which corresponds to the decay of a physical kaon to two pions 
{p = 207 MeV) the value of the derivative of the phase shift with respect to the momentum 
obtained from the phenomenological curve is 9.53 x 10“^MeV~^. Converting this to 
gives —0.216 for the 48^ and —0.221 for the 64^ ensembles. While this makes a signihcant 
difference to the derivative of the phase shift, it represents a relatively small uncertainty 
in the Lellouch-Liischer factor F oc This sum is dominated by the term and 

thus the difference in the Lellouch-Liischer factor between both approaches to calculating 

amounts to 1.1% and 0.6% on the 48^ and 64^ ensembles respectively. 

When quoting our central value we include the Lellouch-Liischer factor evaluated as 
described in the preceding paragraph. In order to estimate the size of the remaining expo¬ 
nential hnite-volume effects we use chiral perturbation theory and include the corresponding 
effects in our systematic uncertainty. Since we are only calculating an estimate, we do not 
use partially twisted chiral perturbation theory, but take both the sea and valence d-quarks 
to satisfy antiperiodic boundary conditions. 

In SU(3 )lxSU(3)r chiral perturbation theory, the leading order (LO) and leading loga¬ 
rithmic next-to-leading order (log) contributions to the (27,1) and (8,8) matrix elements are 
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FIG. 8 : Sample loop diagrams which contribute to finite-volume corrections of (27,1) and 
( 8 , 8 ) K —)■ TTTT matrix elements in chiral perturbation theory. 


given by [ 6 , 24] 


4ia27 1 


4*0; 


27 


fKfl 


{m\ - ml) , 


~2^K ( 1 


m 


K 


m. 

2 2 


/3(m^,m^,m ) + m^ - 


P 

) I3{ml,ml^,ml) + {m% - Smlm'j^ + 2ml) 


5 mj^ 
4 ml 

TT 

4> 


y V v\ I ^ nr/7 ^2 ^2 

X /l(m^,m^,m^) + ( + -m^ 


{ml) + 


-m 


K 


mi 


-Am\ + Ami) Pm\) + ( + Unr^ ) ^{ml) 




4*0, 


88 






4*0, 


88 


/x/2 P 


5 rn 


K 


A mi 


2m\ f3{ml,ml^,ml) + (m^ - 2 * 7 *^) 


1 Tfl ^ 

x/l(m^, * 7 * 2 , ml) + -—f/3{ml, mj^, tt*^) - ( 4 + 


4 *7*^ 


1 *7*|. 
1 ' 2 mi 


i{m 


K) 


l4-8Gmr 

4 *7*2 / ^ ^ 


^'pRn/ 2^ 
-—*7*2 

4 *7*2 ' 


(43) 


(44) 

(45) 


(46) 


At this order * 7 *^ is given by the Gell-Mann-Okubo relation: 3 * 7*2 _ 4 ^ 2 ^ _ 

The functions £(***2) and f3{q^,ml,ml) correspond to diagrams with one and two pseudo 
Goldstone boson propagators respectively as illustrated in Fig. 8 and they are the only 
sources of finite-volume corrections. They are given by (in Minkowski spacetime) 
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l{m^) = 


/3{q,mi,m2) = 


dk^ 


tl 271 k‘^ — w? + ie 

f dkf' i 

^ J 271 (A;2 — mf) {{q + kY — m|) 


OJl + 0^2 


2uJiUJ2{ql - (wi + UJ 2 Y) 


(47) 

(48) 


where the symbol ^ denotes the summation over k in hnite volume or the integration in 
inhnite volume. Ui = -y/^'^Tr^ and a ;2 = \J{q — kY + rnl- The difference between the 
sum and the integral can be calculated using the Poisson summation formula: 


1 

1 ? 




I 


d^k 

(27r)3 


fik) + 

n^O 


/ 


d^k 

(27r)3 


f{k)e 


iLk-n 


(49) 


where the summation on the left-hand side is over all k = ^n, where n is a vector of 
integers. If / is a function that has no singularities on the real axis, then the second term 
on the right-hand side gives the exponential hnite-volume corrections which we are trying 
to evaluate. 


A. Corrections to £{m^) 


With periodic boundary conditions, applying the Poisson summation formula (49) to i, 
writing k in spherical polar coordinates and integrating over the angles, we obtain for the 
difference between the hnite- and inhnite-volume values of £{rrY) [25] 


Ai{m, L) 


2 

m 


5i{mL) 


m Ki{\n\mL) 




E 

n^O 


n 


(50) 


where Ki is a modihed Bessel function of the second kind, n is an vector of integers and the 
sum is over all n Y (0) 0) 0) ^ ^■ 

Since our choice of boundary conditions breaks the isospin symmetry Eq. (50) does not 
give the correct hnite-volume corrections for all the instances of I which appear in Eqs. (44) 
and (46). Specihcally, 7 r°, and 77 satisfy periodic boundary conditions (so that the 
corresponding hnite-volume corrections are indeed given by Eq. (50)) whereas and 
satisfy antiperiodic boundary conditions for which the hnite-volume corrections to i are 
diherent. In the antiperiodic case, we replace f{k) in Eq. (49) by f{k + q), where q = 
(^)(1,1,1). Shifting the integration variable from k to k + q, we hnd that 5i{mL) in Eq. (50) 
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is now replaced by 


5^{mL) 


4 

mL Ini 

n^O 


(51) 


where the index A denotes that the correction is evaluated for a volume with antiperiodic 
boundary conditions in all spatial directions. The difference from the periodic case is the 
additional factor (^—in the summands. The known formulas in Eqs. (44) and (46) 
do not differentiate between different isospin components, and therefore do not specify which 
linear combination of periodic and antiperiodic corrections should be used. Since we are only 
using these formulas for an approximate estimate of the size of the error, we choose to be 
conservative and to include the larger corrections which are those obtained with the periodic 
boundary conditions given in Eqs. (50). The numerical results are presented in Table V and 
as expected the leading contributions come from the loops with a pion propagator. 


B. Corrections to and 


We now consider the contributions from loops with two meson propagators and which are 
proportional to the function /3. We start by discussing the corrections to uik, and 

/9(m^, rriK, for which in Minkowski space the external energy is below the corresponding 
two-particle cut; e.g. in the external energy in the center-of-mass frame is 

ttItt which is clearly smaller than rriK + In such situations the hnite-volume corrections 
are exponentially small. We postpone the discussion of the contribution which does contain 
the two-particle cut, that proportional to /9(mi^, m^), until the following subsection. 

The corrections to /9(m^, mx, ^^ 7 ^) and are proportional to 



d?k €^^''^{001+002) 
(271^) 2(jJi(jj2{q^ - (cui -F 0 ) 2 )^) 


(52) 


with 


Wi = 


m\ 


and 


UJr, = 


q + k 


+ m2 


(53) 


Because of the angular dependence inside the integrals, we evaluate the integrals numerically. 
With the boundary conditions which we are using the corrections with a K~^ and 7r~ are 
equal and opposite to those with the neutral mesons. In the estimate of the uncertainty we 
conservatively do not exploit the cancellation but take the absolute value in each case. 
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We note that care must be taken when using Eqs. (71) and (73) for the hnite-volume 
corrections to {3 in Sec. VIII of [6]. In Eq. (52) above, the two terms in the factor in the 
denominator of the integrand — (cai + a; 2 )^ come with opposite signs. How this arises 
in hnite-volume Euclidean correlation functions is explained in the Appendix following [5]. 
The corresponding terms in the denominator of Eq. (73) in [6] appear (incorrectly) with the 
same sign. 


C. P{mK,m„,m^) 

Kinematically this case is simpler than the two /9 integrals which were evaluated in Sec. V B 
since the external particle {K) is now at rest which eliminates the angular dependence 
from the integral. Furthermore, both internal propagators satisfy antiperiodic boundary 
conditions. In this case however, the integral for /3 has a pole at a;,r = so the Poisson 

summation formula will give both the exponential and powerlike corrections. The power 
corrections are included as the Lellouch-Liischer factor F in Eq. (40) and we do not include 
these in the estimate of the hnite-volume uncertainty. The evaluation of the remaining 
exponential corrections following the approach of [26] is explained in the Appendix. 


D. Combining the finite-volume corrections 

To one-loop order we write the systematic error associated with the hnite-volume correc¬ 
tions in terms of the ratios AAdiog/AdLO- These are given by: 




+m 


P{m\ - ml) [ 12 


■—( 1 


m 


K 


mt 




2 ]^rn\^2ml] AI3{ml,m\,ml) + 


K 


4 


(m^ - + 2ml)A/3{m%, ml, ml)+[- 


1 m% 1 2 1 5 

, 4 < - li’"* + s’"' 


ml ^ ^ ^ \ Ami 4 


4m^ + 4m^ A£(m^) + 7 —^ A£(m^) 


Ai{m 


( 54 ) 
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Quantity 

48^ lattice 

64^ lattice 

L 

5.48 fm 

5.36 fm 

A/(m2) 

14.32 MeV^ 

16.39 Me 

A/(m4) 

(9.05 X 10-4 MeV^ 

(1.03 X 10-3) MeV^ 

A/(m4 

(1.32 X 10-Q MeV^ 

(1.52 X 10-4 MeV^ 

A/3{mn,mK,m^) 

3.0 X 10-^ 

3.0 X 10-^ 

A/3{m^,mK,mP 

5.0 X 10-5 

5.2 X 10“^ 

AP{mK,m^,mP 

6.67 X 10"^ 

6.97 X 10-3 

^^'^(27,1) 

^( 27 , 1 ) 

0.022 

0.024 

Vs,8) 

0.024 

0.026 


TABLE V: Contributions to our estimate of the exponentially suppressed finite-volume 

errors. 


and 


AAff4 


1 

■ P 
1 


^ 7T7^ 

ITTy “ ) APirn'i, mp + (m^ - 2m^)A/3(m4, ml, ml) 

^ ' ' ^'TT 


4 m'i 


A(5{ml, mj^, m. 


4+l!!4) A£(m2 


\ 2 m2 

5 mi- Ox 3 m2 






4 m: 


y 


4 m2 


(55) 


The numerical values of these ratios for the 48^ and 64^ ensembles are shown in Table V. 


VI. THE ERROR BUDGET 

In this section we discuss the two remaining systematic errors: those which arise because 
the meson masses and the two-pion energy are not quite physical and those introduced by 
the perturbative Wilson coefficients. Finally all of the systematic errors in our results for 
the real and imaginary parts of A 2 are summarized in Tables IX and X, respectively. 

The volume, boundary conditions and quark masses have been chosen to enable simula¬ 
tions of physical K —)■ tttt decays. Nevertheless, since the volume and quark masses have to 
be chosen a priori, the output values of the meson masses and two-pion energies will be a 
little different from the physical values (see Table I). In order to estimate the corresponding 
uncertainty we follow the procedure described in [2, 27] and outlined below. We use mea¬ 
surements on 60 quenched conhgurations on a 24^ lattice with a~^ = 1.31 GeV performed 




with three values of the light-quark masses, hve strange-quark masses and the application 
of antiperiodic boundary conditions in ntw = 0, 1, 2 and 3 directions. These measurements 
are used to determine the coefficients in the following phenomenological formulas: 

mly = Bo{mx + rriy) + Bi, (56) 

= AQ{nt^)mi + Ai(ntw), (57) 

A 2 = Co(ntw)ms + Ci(ntw)mi + C2(ntw), (58) 

where mi and mg are the masses of the light and strange quarks, mxy is the mass of the 
meson consisting of x and y valence quarks (which can be either light or strange) and riu, 
is the number of directions in which the antiperiodic boundary conditions would have to 
be imposed on the quenched lattice to get the correct two-pion energy. Note that ntw does 
not have to be an integer, and is given instead by = ntw7r^/T^, where p is the center-of- 
mass momentum of each pion. The full list of coefficients A, B and C obtained from these 
quenched conhgurations was presented in [27] and is reproduced in Table VI. 

We can use the coefficients in Table VI to determine A 2 on the quenched ensembles for 
any choice of {m.^, mx, 77,r7r}- We exploit this possibility for three sets of parameters: (i) the 
physical masses mx = = 493.7MeV, = 139.6 MeV; (ii) the values from the 48^ 

simulation given in the third row of Table I and (hi) the values from the 64^ simulation 
given in the fourth row of Table I. We denote the corresponding three estimates of A 2 by 
^q.phys^ ^q,48 respectively, where the superscript q reminds us that the results were 

obtained on the quenched ensembles. We use the differences and — 

as estimates of the systematic error due to unphysical kinematics. 

The results are: 


Re(A^’P^y") 

= 2.25 X 10-®GeV, 

Im(A^’P^y^) 

= -1.344 X 10-^2 GeV, 

(59) 

Re(Af«) 

= 2.29 X 10-® GeV, 

Im(Af«) 

= -1.341 X 10-^2 GeV, 

(60) 

Re(Af^) 

= 2.36 X 10-® GeV, 

Im(Af^) 

= -1.329 X IQ-^^GeV. 

(61) 


The differences in Eqs. (59)- (61) translate to an estimated 1.8% error on Re(A 2 ) and 0.2% 
error on Im(A 2 ) on the 48^ ensemble and a 4.5% difference for Re(A 2 ) and 1.1% difference 
for Im(A 2 ) on the 64^ ensemble. These numbers are obtained from the difference of the 
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ntw 

0 

1 

2 

3 

^0 

Bo 

Bi 

ReC'o(GeV) 

ReGiioeV) 

ReCajoeV) 

ImGo(GeV) 

ImGi(GeV) 

ImGa^GeV) 

17.53(16) 

0.0273(12) 

2.124 

0.00692 

1.016(55) X 10-^ 
1.697(89) X 10-® 
2.53(51) X 10-9 
-1.06(31) X 10-12 
5.54(79) X 10-11 
-1.689(64) X 10-12 

17.14(73) 

0.1038(60) 

2.124 

0.00692 

1.43(11) X 10-^ 
1.29(18) X 10-^^ 
1.08(12) X 10-® 
-4.6(3.3) X 10-13 
3.39(91) X 10-11 
-1.392(66) X 10-12 

14.9(2.3) 

0.202(18) 

2.124 

0.00692 

1.53(25) X 10-^ 
1.45(38) X 10-6 
1.68(25) X 10-3 
4.4(7.4) X 10-13 
2.1(1.6) X 10-11 
-1.24(12) X 10-12 

24.5(9.5) 

0.196(82) 

2.124 

0.00692 

1.78(54) X 10-^ 
4.22(97) X 10-6 
-2(67) X 10-1° 
2(11) X 10-13 
-1.8(3.2) X 10-11 
-7.5(1.9) X 10-13 


TABLE VI; Parameters used for extrapolations on the 24^ quenched ensembles. 



(27,1) 

(8,8) 

00 

Y 

X 

^LO 

^NLO 

..NLO 

ai 

0.26696 

-0.0035185 

0.290342 

-0.00397252 

4.260055 X 10-6 
-2.026445 X lO-^ 
4.70099 X 10-6 
-8.09555 X 10-6 

-1.0063 X 10-6 
2.447741 X 10-^ 
-5.22390 X 10-6 
3.26016 X 10-^ 


TABLE VII; Wilson coefficients at 3 GeV in the MS scheme at leading order (LO) and 

next-to-leading order (NLO). 


simulated results from those at the physical point (normalized by the result at the physical 
point). These uncertainties are included in Tables IX and X under the label “unphysical 
kinematics.” 

To estimate the error in the Wilson coefficients, we compare the results for A 2 using 
Wilson coefficients calculated at leading order and next-to-leading order. We have used the 
set of coefficients evaluated in the MS scheme at 3 GeV, which are shown in Table VII [28], 
and the standard parametrization of Wilson coefficients was used, i.e. Ci = Zi -h ryi where 
r is the ratio of GKM matrix coefficients r = — Throughout this paper we use the 



LO 

NLO 

Re(A2) 483 
Im(A2) 48° 
Re(A2) 643 

Im(A2) 64° 

1.293(11) X 10-3 
-5.551(45) X 10-13 
1.3410(89) X 10-3 
-6.037(71) X 10-13 

1.386(12) X 10-3 
-6.174(49) X 10-13 
1.4386(95) X 10-3 
-6.548(78) X 10-13 


TABLE VIII; Gomparison of matrix elements calculated with leading order (LO) and 
next-to-leading order (NLO) Wilson coefficients. 
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ReA2 systematic errors 

48^ 

643 

cont. 

NPR (nonperturbative) 

0.1% 

0.1% 

0.1% 

NPR (perturbative) 

2.9% 

2.5% 

2.9% 

Finite-volume corrections 

2.2% 

2.4% 

2.4% 

Unphysical kinematics 

1.8% 

4.5% 

4.5% 

Wilson coefficients 

6.8% 

6.8% 

6.8% 

Derivative of the phase shift 

1.1% 

0.6% 

1.1% 

Total 

8% 

9% 

9% 


TABLE IX: Systematic error breakdown for Re A 2 


ImA2 systematic errors 

483 

643 

cont 

NPR (nonperturbative) 

0.1% 

0.1% 

0.1% 

NPR (perturbative) 

7.0% 

6.2% 

7.0% 

Finite-volume corrections 

2.4% 

2.6% 

2.6% 

Unphysical kinematics 

0.2% 

1.1% 

1.1% 

Wilson coefficients 

10% 

8% 

10% 

Derivative of the phase shift 

1.1% 

0.6% 

1.1% 

Total 

12% 

10% 

12% 


TABLE X: Systematic error breakdown for Im A 2 


particle data group convention for the matrix elements, where Vus = 0.97425, Vud = 0.2252 
and r = 0.0014148 — 0.0005558L The results for matrix elements calculated at leading and 
next-to-leading orders are shown in Table VIII. From the differences between the entries in 
the columns marked as LO and NLO we estimate that the uncertainties are 6.8% for Re(A2) 
on both sets of ensembles and 10% (8%) for Im(A2) on the 48^ (64^) ensembles. 

Tables IX and X show our estimates of systematic errors associated with the results 
for Re(A2) and Im(A2) presented in this paper. The evaluation of the continuum limit of 
A 2 is discussed in the following section. As will be seen, the systematic error associated 
with this extrapolation is negligible with respect to the statistical errors. Consequently 
no discretization error is shown in Tables IX and X. The values in the column marked 
“Cont.” are the errors assigned to our continuum-extrapolated results, and are simply the 
larger of the corresponding entries from the 48^ and 64^ columns. We can see that the 
dominant contribution to the systematic error for both real and imaginary parts of A 2 on 
both ensembles comes from the uncertainty in Wilson coefficients. 
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481 

641 

Phys. 

m^/mn 

0.08296(17) 

0.08220(19) 

0.08073 

mx/mn 

0.29740(32) 

0.29982(37) 

0.29643 


TABLE XI: The ratios of the pion and kaon mass to the Omega baryon mass on the 48^ 
and 64^ ensembles as well as the physical value. 

VII. CONTINUUM EXTRAPOLATION 

In this section we discuss the extrapolation of the results obtained on the 48^ and 64^ 
ensembles to the continuum limit. We divide this discussion into two parts. In the first 
we present the complete physical results for the complex amplitude A 2 in the continuum 
limit. As we will observe, the dominant error in our result comes from the perturbative 
error assigned to the Wilson coefficients. This may be reduced in the future if higher order 
perturbation theory results become available or if lattice step-scaling methods are used to 
allow present perturbative results to be applied at a higher energy scale. Therefore, in the 
second part we determine the continuum limit of the individual matrix elements themselves, 
normalized in the regularization-independent (^, and ( 7 , 7 ) schemes. 

A. Continuum limit of Re(A2) and Im(A2) 

As already mentioned in Sec. VI the quark masses used in these ensembles are very slightly 
larger than their physical values. This is illustrated in Table XI, in which we compare the 
physical and simulated values of the dimensionless quantities mT,/mQ and rriK/rnn, which 
are highly sensitive to the light- and heavy-quark masses respectively. In order to determine 
the values of the lattice spacing we must therefore perform a short chiral extrapolation; this 
is achieved using a simultaneous chiral and continuum “global £t” that incorporates data 
from both ensembles. Since the (renormalized) quark masses on the two ensembles are very 
similar, we must include additional ensembles in order to have a sufficient spread of masses 
for the determination of the chiral dependence. The full set of ensembles and details of this 
procedure can be found in [7]. 

The determination of A 2 presented here was performed using 76 configurations of the 
48^ ensemble, whereas the lattice spacings in [7] were computed using 80. In order to 
preserve the full correlations between the jackknife samples of A 2 and the corresponding 
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superjackknife samples of the lattice spacing, we repeated the global fit analysis using the 
same 76 conhgurations. The details of the binning are also different. In [7] we binned the 
48^ data over 5 successive measurements (100 MD time units) in order to take into account 
the observed autocorrelations in the data, whereas in the present calculation, as explained in 
Sec. II, we construct 19 bins each of 4 conhgurations. These differences lead to determined 
values of the lattice spacings in Eq. (62) below which are a little different from those in [7]. 
For the 64^ ensembles we use the same set of 40 conhgurations for the evaluation of A 2 and 
the same binning as in the global ht in [7]. 

In order to estimate the systematic errors due to the chiral extrapolation and hnite vol¬ 
ume in the determination of the lattice spacings, we have performed our hts using three 
diherent chiral Ansatze: NLO SU(2) chiral perturbation theory, with and without hnite- 
volume corrections (referred to as the ChPTFV and ChPT forms respectively), and a linear 
Ansatz (referred to as the “analytic” form). In practice we found the lattice spacings ob¬ 
tained from all three Ansatze to be consistent to within a fraction of the statistical error 
due to the dominance of the near-physical data, hence we treat these systematic errors as 
negligible. The hnal results for the values of the lattice spacing are 

064 ^ = 2.3584(70) GeV and = 1.7280(41) GeV, (62) 

where the errors are statistical only. 

The lattice matrix elements M* scale as and so small differences in the lattice spacing 
become amplihed. We have performed the continuum extrapolation of A 2 using the lattice 
spacings obtained with each of the three chiral Ansatze; the extrapolated values are given 
in Table XII. In Fig. 9 we show the continuum extrapolation in the (^, ^) scheme using the 
lattice spacings obtained with the GhPTFV chiral Ansatz. We use results obtained with 
this Ansatz as our central values for each lattice spacing and for the extrapolated value in 
the continuum. 

We obtain an estimate of the component of the chiral extrapolation error arising from 
the lattice spacing determination by taking the difference between the continuum values 
obtained using the GhPTFV and analytic lattice spacings. The full jackknife differences 
are 0.3(2.6) x 10“^° and 0.1(1.2) x 10“^^ for the real and imaginary parts respectively. As 
with the lattice spacings, we cannot resolve these differences within the statistical error; 
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FIG. 9: The continuum extrapolation of Re(y42) (left) and Im(742) (right). The points at 
hnite lattice spacing are taken from Tab. Ill for the (^, ^) intermediate renormalization 

scheme. 


Ansatz 

Re(A 2 ) (xlO-® GeV) Im(A 2 ) (xlO'^^ GeV) 

GhPTFV 

GhPT 

analytic 

1.501(39) -6.99(20) 

1.494(38) -6.96(19) 

1.494(43) -6.96(21) 


TABLE XII; The continuum values of Re(A 2 ) and Im(A 2 ) determined using the lattice 
spacings obtained with each of the three chiral Ansatze. 

hence we set the chiral error to zero. On the other hand the jackknife differences between 
the ChPTFV and ChPT Ansatze are resolvable as they differ only in small Bessel function 
corrections and are thus highly correlated: we obtain 3.4(1.0) x 10“^^ and 1.59(47) x 10”^^ 
for the real and imaginary parts respectively. Nevertheless, these errors are only 5%-8% of 
the statistical error and can therefore also be neglected. This leads to the result 

Re(A 2 ) = 1.501(39) X 10“® GeV and Im(A 2 ) =-6.99(20) x 10"^^ GeV, (63) 

where the errors are statistical. 

Our hnal result for A 2 is obtained by assigning the 9% and 12% systematic errors from 
Tables IX and X as the systematic errors to be associated with the values for Re(A 2 ) and 
Im(A 2 ) given in Eq. (63): 


Re(.42) = 1.60(4).„(14),,.. x lO"® GeV; Im(Xl 2 ) = -6.99(20)...i(84).,.. x lO"'® GeV. 


(64) 
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Systematic errors in ImA 2 /ReA 2 

48^^ 

643 

cont 

NPR (nonperturbative) 

0.1% 

0.1% 

0.1% 

NPR (perturbative) 

7.6 % 

6.7 % 

7.6 % 

Finite-volume corrections 

3.5 % 

3.5 % 

3.5 % 

Unphysical kinematics 

1.8 % 

4.6% 

4.6% 

Wilson coefficients 

12.0 % 

10.5 % 

12.0% 

Derivative of the phase shift 

0 

0 

0 

Total 

14.7% 

13.7% 

15.3% 


TABLE XIII; Systematic error breakdown for ImA 2 /ReA 2 . 


In order to estimate the unknown quantity ImAo, we combine our results for A 2 with the 
experimental values of ReAo = 3.3201(18) x 10“'^ GeV and e'/e = (1.65 ± 0.26) x 10“^ [29]. 
To this end we start by evaluating the ratio ImA 2 /Re^ 2 ; taking into account any statistical 
correlations between the real and imaginary parts by performing the analysis within the 
jackknife procedure. On the two ensembles we hnd 

=-4.45(5)stat(65)syst X 10-'’ and =-4.55(5),tat(62)syst x 10“^ 

V Jve/i2 /4g3 \ K,e/i2 /g43 

(65) 

The systematic errors for this ratio are given in Table XIII; they are generally combined 
in quadrature except for that due to the derivative of the phase shift because the Lellouch- 
Liischer factor cancels in the ratio. It is interesting to note that if instead of adding the 
errors in the Wilson coefficients for ReA 2 and ImA 2 in quadrature as in Table XIII, we had 
calculated the ratios with the coefficients at leading and next-to-leading order respectively 
and taken the difference as a measure of the uncertainty we would have obtained a much 
smaller answer (3.6% instead of 12%). Since the operators which give the dominant contri¬ 
butions to the real and imaginary parts are different, and in the absence of an understanding 
which might suggest a correlation between their Wilson coefficients, we prefer to be cautious 
and take the larger uncertainty. We hnd a similar feature in the NPR perturbative error. 

The continuum extrapolation of the dimensionless ratio ReA 2 /IuiA 2 is milder than that 
of ReA 2 and ImA 2 separately and we obtain 


/ ImA2 A 
\ ReA2 / f 

\ ^ / continuum 


-4.67(72) X 10"^ 


( 66 ) 
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FIG. 10: Dominant contractions contributing to Re(y42): Ci (left) and C 2 (right). 
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FIG. 11: Cancellation of dominant contributions to Re(R 2 ) on the 48^ ensembles with a 
K-tth separation of 27 and the 64^ ensembles with separation 36. 


Using this ratio, we can calculate the electroweak penguin contribution to e'/e, given by 


EWP 


u lmA2 

V2 lei ReR 2 


= -6.6(10) X 10-^ 


(67) 


where we have used the values cn = = 0.04454(12) and |e| = 2.228(11) x lO"^ from [2]. 

This value for (e'/e)EWP is consistent with our previously quoted value —6.25(44)(119) x 
10“^ [2]. Finally, for ImRo we hnd 


ImRn = ReAn 


ImA2 \/2| 


e e 


ReAo 


ui e 


= -5.40(64) X lO-'^GeV 


( 68 ) 


The results in Eqs. (67) and (68) were obtained using our result for ImA 2 /ReA 2 in 
Eq. (66). If instead we take ImA 2 from our calculation, Eq. (64), and combine it with the ex¬ 
perimental result ReA 2 = 1.4787(31) x 10“^ GeV we obtain, ImA 2 /ReA 2 = —4.73(58) x 10“^, 
(eVe)EWP = -6.69(82) x 10"^ and Im Aq = -5.42(63) x lO’^^ GeV. 
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B. Continuum limit of the RI-SMOM matrix elements 

From the error budget in Table XIII we see that the dominant uncertainty is due to 
the Wilson coefficients, which we take to be the difference between the leading and next- 
to-leading order contributions as dehned in [30], where the calculations were based on [31- 
33]. In case the Wilson coefficients in the RI-SMOM schemes become known with better 
precision in the future, we present in Table XIV the 7i~^7r^ matrix elements 

dehned in Eq. (7), with the operators Qi in Eqs. (8) - (10) renormalized in the (^, and (7,7) 
renormalization schemes at a renormalization scale of 3 GeV. These matrix elements together 
with the new Wilson coefficients would enable an improved evaluation of A 2 , without the 
need to recompute the matrix elements. The systematic errors for the (27,1) operator are 
estimated using the entries in Table IX with the NPR(perturbative) and Wilson coefficient 
errors set to zero. This gives the errors of 2.8%, 5.1% and 5.2% for the 48^ and 64^ ensembles 
and in the continnum limit respectively. For the (8,8) operators using the entries in Table X, 
the same procedure leads to systematic errors of 2.6%, 2.9% and 3.0% for the 48^ and 64^ 
ensembles and in the continnum respectively. 

For completeness we also convert these three matrix elements into those 

in the original 10 operator basis as dehned in [34]: 

= 3Mf ^ = 3Mf ^ = 2Mf ^ = 2]%/^^ (69) 

Af,?;, = and A^^., = 2Mf (70) 

where = ((7r7r)/=2 I Qi \ K^)- 

VIII. CONCLUSIONS 

Before briehy snmmarizing our results and discussing prospects for fnture calculations 
we conhrm our Ending, Erst presented in [3], that there is a significant cancellation between 
the two dominant contributions to Rea42- As explained above, Re(a42) is dominated by the 
matrix element of the (27,1) operator and is proportional to the snm of the two contrac¬ 
tions Cl and (7*2 in Fig. 10. While naive factorization, freqnently used for phenomenological 
estimates, suggests that Gi = 3 G 2 because of the color suppression in C 2 , we find a strong 
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Ensemble 

Scheme 

Vki) (GeV*) 

(G'V=>) 

(G'V=>) 

483 

Hi) 

0.04761(39)(133) 

0.7026(52)(183) 

3.892(28)(101) 

643 

Hi) 

0.04848(32) (247) 

0.8412(88)(244) 

4.140(44)(120) 

483 

(7,7) 

0.04473(37)(128) 

0.7112(53)(185) 

3.471(26)(90) 

643 

(7,7) 

0.04664(31)(238) 

0.8477(88) (246) 

3.724(40)(108) 

Continuum 

Hi) 

0.0506(13)(26) 

1.003(22)(30) 

4.43(12)(13) 

Continuum 

(7,7) 

0.0489(13)(25) 

1.007(23)(30) 

4.02(10)(12) 


TABLE XIV: Results for the —)■ { 7171 ) 1=2 matrix elements (defined in Eq. (7)) in 

two non-exceptional RI-SMOM renormalization schemes at the scale 3 GeV. The first error 
is statistical, while the second one is the systematic uncertainty estimated as described in 

the text. 

cancellation between these two contributions. For the 48^ and 64^ ensembles studied in this 
paper, we illustrate this cancellation in Fig. 11. (In Sec. Ill we explain that the numerical 
results in this paper were obtained from correlation functions with even values of The 
choice of = 27 for the 48 ensembles in Fig. 11 is made to ensure that the cancellation 
is illustrated at the same value of t-j^n m physical units on the two sets of ensembles.) As 
explained in [3] we believe that this cancellation is a signihcant component in explaining the 
AI = 1/2 rule. Although we have not completed the calculation of Aq at this stage, we note 
that the contributions of the (27,1) operator all contribute with the same sign. A similar 
partial cancellation occurs between the two corresponding contractions in the evaluation of 
the Bx parameter of neutral kaon mixing as pointed out in [35] and subsequently conhrmed 
in [3, 36]. 

Our ab initio determination of A 2 shows clearly that phenomenological approaches based 
on the dominance of naive factorization are not consistent. We note however, that there 
were nonlattice studies based on chiral perturbation theory and the 1/N expansion, where N 
is the number of colors, which indicated that C 2 may have the opposite sign to Ci [37, 38]. 
Of course, as illustrated in our results above, the 1/N expansion per se is not a good 
approximation; C 2 is suppressed by 1/N and yet is comparable to Ci. In different ways, the 
authors of [37, 38] combine the expansion with leading short- and long-distance logarithms. 
In [37] the authors use an Ansatz for matching the perturbative short-distance contributions 
and long-distance effects based on a chiral Lagrangian for mesons. In [38] the authors 
compare the experimental value of ReA 2 with the leading term of the expansion to deduce 




38 


that C 2 should be negative. For recent discussions of these two early approaches, stimulated 
by our lattice QCD result [1, 2] and written by subsets of their original authors, we refer 
the reader to [39, 40]. 

Our earlier calculation of A 2 was performed on an ensemble at a single coarse lattice 
spacing, a~^ = 1.364 GeV [1, 2], and so not surprisingly the dominant systematic uncertainty 
was due to discretization errors. We estimated these to be about 15%, although with only 
a single lattice spacing this could only be an estimate. In the present paper we repeat and 
rehne the earlier calculation which is now performed on two hner ensembles with different 
lattice spacings, allowing for a continuum extrapolation. We have determined Re %2 to be 
1.50(4)stat(14)sys X 10“® GeV. This is consistent with the experimental values of 1.4787(31) x 
10“® GeV from charged kaon decays and 1.570(53) x 10“® GeV from neutral kaon decays. We 
have also calculated the imaginary part of A 2 to be —6.93(20)stat(84)sys x 10“^^ GeV, which 
was unknown until [1, 2]. [We recall that the corresponding results from our earlier work 
were Re2l2 = l-38(5)stat(26)sys x 10“® GeV and ImR2 = — 6.54(46)stat(120)sys x 10“®GeV.] 
Our results for Im and Re A 2 imply {e'/e)EWP = —6.6(10) x 10“"^. This can be compared 
to the result obtained via hnite energy sum rules [41], Re(e'/e)£;wp = —(11.0 ± 3.6) x 10“^ 
(see also results based on vacuum saturation [41, 42]). We also mention for completeness 
that the continuum value of the two-pion phase shift is <5 = —0.203(43).’ 

The errors are currently dominated by systematic uncertainties, the largest of which is due 
to the uncertainty in the (perturbative) evaluation of the Wilson coefficients (see Tables IX 
and X). It is testimony to the huge progress in the precision of lattice calculations that this 
is the case. We have aimed to be conservative in estimating this error, taking the difference 
between the lowest order and the next-to-lowest order as the uncertainty. The natural way 
to decrease this error is to perform higher-order perturbative calculations in the standard 
model but it may also be possible to use step scaling to increase the renormalization scale 
in the intermediate schemes (such as the RI-SMOM schemes used in this study) and hence 
to increase the scale at which the matching to the MS scheme is performed and at which 
the Wilson coefficients are calculated. It will be interesting to explore this possibility. 

In order to have a fully quantitative understanding of the AJ = 1/2 rule, to determine 
e'/e and to compare the result to the experimental value e'/e = (1.65 ± 0.26) x 10“^ we 
need to perform the evaluation of Aq at physical kinematics. A key ingredient which makes 
the calculation of A 2 feasible is the use of the Wigner-Eckart theorem described in Sec. III. 
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Together with the choice of volume and the use of antiperiodic boundary conditions for the 
d-quark in all three spatial directions, it ensures that the energy of the two-pion gronnd state 
is equal to uik- Unfortunately this approach cannot be directly applied to the calculation 
of Af)] in particular the breaking of isospin symmetry by the bonndary conditions invali¬ 
dates the calculation. For example, the 7r° remains at rest with the antiperiodic boundary 
conditions, whereas the charged pions have nonzero momentum. More sophisticated bound¬ 
ary conditions mixing qnarks and antiqnarks and an isospin rotation, the so-called G-parity 
bonndary conditions [11, 43-46], must therefore be used instead for both the valence and 
the sea qnarks. The evalnation of Aq with G-parity bonndary conditions is well nnderway 
and exciting progress has recently been reported in [4] and we anticipate the hrst complete 
calculation of Aq, albeit on a single lattice spacing, within the next year. 
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FIG. 12: Contribution to the correlation function in which two pions are produced by an 
operator at the origin (grey circle), and rescatter by the strong interactions denoted by the 

hlled circle. 
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Appendix A: P{mK,Tn.jr,Tn„) and the Lellouch-Liischer Factor 


In Sec. V we use chiral perturbation theory to estimate the hnite-volume corrections in our 
calculation of A 2 and consider the differences between the hnite-volume sums and inhnite- 
volume integrals in and P{q,mi,m 2 ) dehned in Eqs. (47) and (48). In the case with 

q = {rriK-i 0) and mi = m 2 = m^r, (3{mK, m.^) in Minkowski space has an imaginary part 

which leads to hnite-volume corrections in Euclidean space which decrease only as inverse 
powers of the volume and not exponentially. These power corrections are the one-loop 
chiral perturbation theory (NLO ChPT) contributions to the Lellouch-Liischer factor F in 
Eqs. (40) and (41). This factor is included fully in our analysis and so we must not include 
it again from NLO GhPT. A detailed study of how the Lellouch-Liischer factor arises in 
one-loop GhPT was performed in [5], but we hope that it will be useful to summarize the 
main points here. 

In Minkowski space, performing the /cq integration in the center-of-mass frame we obtain 




(Fk 1 

(27r)3 _ Auj‘^(k) + is] ’ 


where uj‘^{k) = 


mi 


(Al) 
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In finite-volume Euclidean space we evaluate the correlation function illustrated in Fig. 12. 
The kaon propagator is irrelevant for our discussion and so we amputate it, and consider 
the two pions to be created at the origin, to rescatter and to be annihilated on the time 
slice at ty. After performing the integrals over yi, y 2 (with phase factors and 
respectively) and x and exploiting the resulting 5 functions, we obtain for this contribution 
to the correlation function: 


I = 



dEi 

El + 


^i{Ei-E2)t ^i(E3-Ei)(ty-t) 


(A2) 


where in a finite volume the integral over k is replaced by the corresponding sum. Here 


-|- and oj‘1 = a;| = u‘^{q) = -l- so that are not 


uf = U 2 = = 

integration variables. 

The energy integrals can now be performed by contour integration; there are three con¬ 
tributions depending on the value of t. 


1. The first contribution is from the interval — cxd < f < 0 and gives 

^ ^ ^-2u^{q)ty r ^ 3 ^ ^ 

^ “ 32a;2(^ J (27r)3 ^ 2 ^^) {uj{k) + u{^) ’ 


(A3) 


2. The second contribution comes from the region 0 < t < ty and gives 


h = 


^-2u,(q)ty r ^ 3 ^ 




32a;2(^ J u‘^(k) {u{k) — uj{q)) 


(1-e 


(A4) 


3. Finally we have the contribution from the region ty < t < 00 which gives 


h = 


d^k 


3 — 2 Cj(fc) ty 


32 {2n)^ oo^k) {uj{k) + uj{^) ' 


(AS) 


The contribution to the amplitude is given by the coefficient of 


Q-2oj(q)ty 

4:Uj‘^{q) 


In hnite volume (FV) the integrals over k are replaced by the corresponding sums and we 
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obtain the following three contributions. The first two are 




8L3^ uj2(k) {u{k)+u{^) 


from the region t < 0 , and 




+ 


L3 J 8L3 ’ 


(A 6 ) 


(A7) 


from the region 0 < t < where Ug is the degeneracy of states with k = q. The term 
proportional to ty is the FV correction to the two-pion energy and it can be checked that 
this is correctly given by the Liischer quantization condition [5]. Finally from the region 


ty < t < oo we have 


n= (4 


(A 8 ) 


.LV IQu^q) ■ 

We now separate the terms with \k\ = \q\ from those where \k\ ^ |g |. When \k\ = |g |, 
we hnd a contribution 

z/q 1 f 1 1 ] 

Z3 4:;^\4^^4^J ’ ^ 

where the first term in the braces corresponds to Ti and the second corresponds to T 3 . The 
contribution from T 3 is cancelled by the FV correction to the matrix element of the two-pion 
interpolating operator at ty [5] whereas the one from Ti is a contribution to the FV effects 
in the amplitude. 

The contributions from \k\ 7 ^ |^| come from Ti and T 2 and can be combined to give 


1 


E 


KW -‘^"( 9 )) 


(AlO) 


Thus in Euclidean hnite volume we obtain 




16L3E3 ’ 


(All) 
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where it is convenient to define 




itkl 


(A12) 


and the corresponding integrals by 


f d^k 1 

4 J (2ir)3ijn(fc)(;j2(^_(^2(j5)' 


(A13) 


Relating this sum to the corresponding integral gives the Lellouch-Liischer factor [5]. 
We now make this more specific and determine the exponentially small corrections. In 
the difference S[ — Sq there is no term with a pole at u{k) = u{q) so that this difference 
can be related to the corresponding integral using the Poisson summation formula and the 
exponentially small finite-volume corrections can be identified: 


S[ - S’, = 


E 


u{k){u{k) + u{^) 


(A14) 


1 


E 


AL^u{^ ^ uj{k){uj{k) -h uj{^) 8L^u^{^ 


— Ji — Jo + 


Va 


8LJoj'^{q) 


+ Cl ,0 


(A15) 


Thus we see that the finite-volume and infinite-volume results are related by 




Vr, 


16L3E3 


— Ji — Jo + S'q + 


3z/„ 


16^3^3 


+ ^ 1,0 , 


(A16) 


where represent the exponentially small corrections, 


ei,o 


1 

8'k‘^uj{^L 



sin(n/cL) 

uj{k){ijj{k) + oj{q)) ’ 


(A17) 


and n and k are \n \ and \k \ respectively. It was shown in [5] that + + is precisely 

the one-loop contribution to the Lellouch-Liischer factor. The residual exponentially small 
finite-volume effects are given by (The ultraviolet divergence cancels in the difference 
Jq ~ S'q, but if the zeta function regularization is used, as in [22], then J, = 0.) 

We have presented the above detailed discussion because we believe that there is a misun- 
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derstanding in the literature. In Eqs. (71) and (73) of [6], the authors take the hnite-volume 
corrections in /9(mx, m^) in Euclidean space to be the difference between the momentum 

integral and the corresponding sum over the integrand in Eq. (Al) but with the replacement 
nij^ — 4:u‘^{k) —)■ + 4a;^(/c) in the denominator. Since there would now be no singularity 

in the denominator, the hnite-volume corrections would be exponential and there would be 
no Lellouch-Liischer factor. The above derivation demonstrates instead the origin of the 
power corrections in the volume. 

Throughout the above discussion we assumed periodic boundary conditions in all three 
spatial directions so that ki = rii x [271/L) where n* is an integer. In our determination of 
A 2 we use antiperiodic boundary conditions in all three directions so that 


ei,o 


1 

Sn‘^uj{q)L 


E 


(- 1 ) 


rix+ny+n^ 


n 



sm.{nkL) 

u{k){(jj{k) + oj{q)) 


(AI8) 
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